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We present a complete systematically theoretical study of multifragmentation for asymmetric 
colliding nuclei for heavy- ion reactions in the energy range between 50 MeV/nucleon and 600 
MeV/nucleon by using soft and hard equations of state. This study is performed within an isospin- 
dependent quantum molecular dynamics model. To see the effect of mass asymmetry, simulations 
are carried out in the absence of Coulomb interactions. Coulomb interactions enhances the produc- 
tion of fragments by about 20%. We envision an interesting outcome for large asymmetric colliding 
nuclei. Although nearly symmetric nuclei depict a well known trend for rising and falling with a 
peak around E= 100 MeV/nucleon, this trend, however, is completely missing for large asymmetric 
nuclei. Therefore, experiments are needed to verify this prediction. 

PACS numbers: 25.70.Pq, 25.70.-z, 24.10.Lx 



I. INTRODUCTION 

The ability to understand the properties of nuclear 
matter at the extreme conditions of temperature and 
density is one of the challenges in present-day nuclear- 
physics research. Such extreme conditions can be gener- 
ated in a heavy-ion-induced reaction at intermediate en- 
ergies . The outcome of a reaction depends on the 
incident energy, the impact parameter, as well as on the 
asymmetry of the colliding partners [lHa] . For symmet- 
rically heavy colliding nuclei at central impact parame- 
ters, two primary fragments are formed: one that is the 
projectilelike fragment and the other that is the target- 
like fragment. These excited fragments deexcite through 
various exit channels: evaporation of light particles and 
emission of intermediate- mass fragments (IMFs) [lj-Q. 
The excitation energy deposited in the system at low in- 
cident energies is too small to allow the break up of the 
nuclei into fragments. With an increase in the incident 
energy, colliding nuclei may break into dozens of frag- 
ments consisting of light, medium, and heavy fragments. 
The size of the fragments and physics behind their for- 
mation differs in different physical conditions. No such 
fragments will survive at extremely high incident ener- 



gies. 

Among various theoretical models developed to study 
these reactions, one can group them into those, which 
are statistical in nature [6[, and others, which take 
the dynamics of the reaction into account and, hence, 
are capable of investigating the evolution of the_frag 
mentation and nucleon-nucleon correlations [lj- 
Interestingly, both types of models (although different 
in their assumptions) are able to explain one or the 
other feature of the experimental findings. Here, we 
will only concentrate on the dynamical model. A care- 
ful analysis of experimental efforts reveals either that 
one has studied the collision of symmetric nuclei (e.g., 
rgAu 197 +79 Ait 197 ) [3j or that one has studied asymmet- 
ric colliding nuclei(e.g., 2oAr 40 + 2 iSc 45 , 2oAr 40 +7gAu 
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eC 1 ^ +79 Au w ', 20 Ca 4U + 79 Au Li ") 3]. Although the 
dynamics for symmetrically heavy nuclei is prominently 
exposed in experimental and theoretical studies, little 
attention is paid to the collision of asymmetric nuclei. 
We, at the same time, know that the symmetry and the 
isospin play decisive roles in a reaction. We want to dis- 
cover the fragmentation of asymmetrically colliding pairs 
in the following different ways, (i) In the first case, we 
will perform a systematic study of the emission of various 
fragments as a function of the asymmetry rj of a reaction. 
The asymmetry of a reaction can be defined by the asym- 
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metry parameter i] = (Ap — Ap)/(Ap + Ap) 
and Ap are, respectively, the masses of the target and the 
projectile, rj — corresponds to the symmetric reaction, 
whereas nonzero values of rj define different asymmetries 
of the reaction. It is worth mentioning that the out- 
come and the physical mechanism behind the symmetric 
and asymmetric reactions are entirely different 0, [^, [lol | . 
Here, for systematic analysis, we start from the symmet- 
rically colliding partners (77 = 0), and then, asymmetry 
parameter rj is varied gradually (rj = -0.8 to 0.8) by keep- 
ing the total mass of the system fixed. Such an experi- 
ment was performed by Betts 11] in 1981, where fusion 
probabilities were measured for different colliding pairs, 
which lead to the same compound nucleus, 
(ii) In the second case, the projectile mass is varied 
from 16 to 56 units, while the total mass of the sys- 
tem is kept fixed. For example, we study the reactions 

_r n!6 1 V„136 c,-28 , y .124 c32 , c„120 
Ot $U +54 Ae , i 4 5l +54 Ae , i 6 D +50 Oil , 

20 Ca 40 +50 Sn 112 , 2i Cr 50 + 44 Ru 102 , 26 Fe 5e + 44 Ru 96 , 



etc. The target isotope is chosen to be a stable one. 
Since, we will neglect the Coulomb effect, we may say 
that it leads to the same compound nucleus. 
In a recent communication [12j |. Liu studied the isospin 
effects on the process of multi-fragmentation and dis- 
sipation by considering the two pairs of colliding sys- 
tems Zn 76 + Ar ia and Kr 76 + Ca 40 , Cd 120 + Ar i0 and 
Xe 120 + Co 40 for central collisions. Another study 0] 
focused on the isospin effects of the mean-field and two- 
body collisions on nucleon emissions at intermediate en- 
ergies. This study showed that the neutron-proton ra- 
tio of preequilibrium nucleon emission and the neutron- 
proton differential and elliptical flows are the probes for 
extracting the isospin-dependent mean field at a lower 
beam-energy region. Because of less compression in 
asymmetric reactions, most of the deposited excitation 
energy is in the form of thermal energy. Our present 
study will shed some light on the effect of the asymme- 
try of a reaction on fragmentation, where a great amount 
of energy is in the form of thermal energy. 



The present analysis will be carried out within the frame 
work of the isospin-dependent quantum-molecular dy- 
namics (IQMD) model 0,0. Our paper is organized 
as follows: We briefly discuss the model in Sec. II. Our 
results are given in Sec. Ill, and we summarize the results 
in Sec.IV. 

II. THE MODEL 

The IQMD model treats different charge states of 
nucleons, deltas, and pions explicitly 14j, as shown in 



the Vlasov-Uehling-Uhlenbeck (VUU) model [l5j]. The 
IQMD model was successfully used in analyzing the large 
number of observables from low to relativistic energies 
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13-15]. One of its versions (quantum-molecular 



dynamics) has been very successful in explaining the 



subthreshold particle production 116 
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the collective flow 0, [17J , the disappearance 



the multifragmen- 



tation 

of flow [17j. and the density temperature reached in 
a reaction [8j. We will not take relativistic effects 
into account, since there is no effect 0] in the energy 
domain in which we are interested. The isospin degree 
of freedom enters into the calculations via both cross 
sections and mean field 0]. The details about the 
clastic and inelastic cross sections for proton-proton and 
neutron-neutron collisions can be found in Refs. HQ. 
In this model, baryons are represented by Gaussian- 
shaped density distributions 



fi(r,p,t) 



-('— r-j{t)) 2 -(P-Pj(t)) j .2L 



TT 2 H 2 



(1) 



Nucleons are initialized in a sphere with radius R = 
1.12 A 1 / 3 fm, in accordance with the liquid- drop model. 
Each nucleon occupies a volume of h 3 so that phase 
space is uniformly filled. The initial momenta are 
randomly chosen between and Fermi momentum pp. 
The nucleons of the target and the projectile interact 
via two- and three-body Skyrme forces and the Yukawa 
potential. The isospin degrees of freedom are treated 
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explicitly by employing a symmetry potential and the 
explicit Coulomb forces between protons of the colliding 
target and protons of the projectile. This helps to 
achieve the correct distribution of protons and neutrons 
within the nucleus. 

The hadrons propagate by using Hamiltonian equations 
of motion: 



(2) 



dfl d < H > dpi d < H > 

dt dpi dt dfl 

with 

<i?>=<T> + <y>as the Hamiltonian. 
= E|| + EE / fi(r,P,t)V ij (?,f) 

i i j>i 

fj(P,p',t)dfdPdpdp. (3) 

The baryon-baryon potential V %3 , in the preceding 
relation, reads as 



V lj {r' -f) 



bkyrme ' Yukawa 



v Coul + v Sym 



hS(P -f) + t 2 S(P - f)pi- 1 



exp(\ r' -f\/p) ZjZje 2 
(| r' — r\/fi) | r' — f \ 



ti-TpiWi-f}). 

Po 



(4) 



Where p. = 0.4/m, t 3 = -6.66^/6^, and t 4 = WOMeV. 
The values of t\ and t 2 depends on the values of a, (3, 
and 7 Here, Zj and Zj denote the charges of the 
i th and j th baryons, and Tg, Tg are their respective T3 
components (i.e. 1/2 for protons and -1/2 for neutrons). 
The Meson potential only consists of the Coulomb inter- 
action. The parameters [i and t±, ,t± are adjusted 

to the real part of the nucleonic optical potential. For 
the density dependence of the nucleon optical potential, 
standard Skyrme-type parametrizations are employed. 
The Skyrme energy density has been shown to be 
very successful at low incident energies, where fusion 
is the dominant channel 0, [ltj]. The Yukawa term 
is quite similar to the surface-energy coefficient used 



in the calculations of the nuclear potential for fusion 
[ijj]. The choice of the equations of state (EOS) (or 
compressibility) is still a controversial one. Many studies 
advocate softer matter, whereas, a greater number of 
studies indicate matter to be harder in nature 00. 
We will use both hard (H) and soft (S) EOS that have 
compressibilities of 380 and 200 MeV, respectively. 
The symmetry energy is taken into account by introduc- 
ing 



V» m =U-TiTi.6(r' i -f' j ). 

Po 



(5) 



The binary nucleon-nucleon collisions are included by 
employing the collision term of the well-known VUU- 
Boltzmann-Uehling-Uhlenbeck equation 
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The 



binary collisions are allowed stochastically, in a similar 
manner as performed in all transport models. During 
the propagation, two nucleons are supposed to suffer a 
binary collision if the distance between their centroids 
is, 



< 



Utot 
7T 



vtot = a(y/s,type), (6) 



where "type" denotes the in-going collision partners 
(N-N, N-6, N-7T...). In addition, Pauli blocking (of the 
final state) of baryons is taken into account by checking 
the phase space densities in the final states. The final 
phase space fractions Pi and P 2 which are already 
occupied by other nucleons, are determined for each of 
the scattering baryons. The collision is then blocked 
with probability 



block 



1 - (1 - Px)(l - P 2 ). 



(7) 



The delta decays are checked in an analogous fashion 
with respect to the phase space of the resulting nucleons. 
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III. RESULTS AND DISCUSSIONS 



In the present calculations, a simple spatial clusteri- 
zation algorithm dubbed as the minimum-spanning-tree 
method is used to clusterize the phase space fl]. We, 
however, also acknowledge that more microscopic algo- 
rithmic routines are also available in literature 0, 0] . By 
using the asymmetric (colliding) nuclei, the effect of mass 
asymmetry can be analyzed without varying the total 
mass of the system. As noted previously, the experimen- 
tal studies by the Michigan State University, miniball and 
ALADIN [21] groups vary the asymmetry of the reaction, 
whereas the plastic ball and FOPI experiments are 
only performed for symmetric reactions. 
The effect of mass asymmetry on fragmentation is 
demonstrated in Fig.l. Here, relative multiplicity Rm 
is defined as = | ( - M ^ t ~ s Ms - ) | where M A and M s are the 
multiplicities of various fragments obtained in the asym- 
metric and symmetric colliding nuclei, respectively. The 
relative multiplicity of free nucleons, light mass fragments 
(LMF's) (2 < A < 4), medium mass fragments (MMF's) 
(3 < A < 8), and IMFs (5 < A < A tot /6) follows hyper- 
bolic behavior. Here, the total mass was kept fixed at 
140 units, and the mass of the target and the projectile 
was varied in steps of 10 units (e.g., Ap = 130, At = 10, 
Ap = 120, At — 20 etc). As asymmetry 77 shifts toward 
the positive side, the target fragmentation takes place 
since At >> Ap. In contrast, at rj = -0.8, the projec- 
tile fragmentation takes place since Ap » At- Since 
the emission of the nucleons (protons and neutrons) is 
maximum in the participant zone, one sees lesser nucleon 
emission compared to the emission of lighter fragments 
such as LMFs and MMFs. In this study, all reactions are 
performed in the laboratory frame. Note that the mirror 
reactions are also studied by the FOPI collaboration [3] . 
Now, we confine our study to particular asymmetric 
systems such as 8 16 +54 Xe 136 , 14 Si 2S + 54 Ae 124 , 
16 S 32 +5 Sn 120 , 20 Ca 40 +5 5n 112 , 24 Cr™ + 44 Ru 102 , and 
26-Fe 56 +44 Ru 96 at different incident energies. To see the 
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FIG. 1: (Color online) The relative multiplicity of different 
fragments as a function of mass asymmetry r\ at E = 150 
MeV/nucleon and impact parameter at b = 0.3. 



effect of the Coulomb interactions, in Fig. 2, we display 
the final phase space of a single event of 26^ e56 +44 Ru 96 
(77= 0.2) (upper panel) and 14 Si 28 +54 AV 24 (77=0.6) 
(lower panel) at a fixed center-of-mass energy of 250 
MeV/nucleon with and without the Coulomb interac- 
tion. Here, the phase space of LMFs [(2 < A < 4)], and 
IMFs [(5 < A < At t/6)] is displayed. Irrespective of the 
Coulomb interaction, the reaction with 77 = 0.2 leads to 
isotropic emission compared to the reaction with 77 = 0.6 
that projects a nearly binary character. Since LMFs orig- 
inate from the midrapidity region, they are better suited 
for studying the effect of asymmetry on the reaction dy- 
namics. 

Generally, obvious effects associated with the asymmetry 
of a reaction are caused by thr Coulomb interaction. To 
understand the role of asymmetry beyond the Coulomb 
effects, we switch off the Coulomb force in further analy- 
sis. Additionaly, we keep the center-of-mass energy fixed 
throughout the analysis. 

In Fig. 3, we compare the effect of the Coulomb forces 
on the multiplicities of various fragments at E c , m ~ 50 
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FIG. 2: (Color online) The final-state phase space of a single 
event for the reaction of 26 Fe 56 +44 -R« 96 with and without 
the Coulomb effects. Here, center-of-mass energy is E c . m . — 
250 MeV/nucleon and impact parameter is b = 0.3. Different 
symbols denote LMFs and IMFs. 

MeV/nucleon and E c m ,— 250 MeV/nucleon. The asym- 
metry of the reaction is varied by using projectiles with 
masses between 16 and 56. We see a clear effect of 
Coulomb forces. As expected, it is maximum for the 
low incident energies, and this effect diminishes as we 
move toward higher incident energies. An enhanced ef- 
fect emerges at larger asymmetric reactions. 

The production of the heaviest fragment A max , the 
free nucleons, the LMFs (2 < A < 4), and the IMFs 
(5 < A < A tot /6) shows expected behavior. The heavier 
mass continues to grow, and it is close to the mass of the 
reacting partners for larger asymmetries. In contrast, 
the production of free nucleons, LMFs, and IMFs shows 
a reverse trend with the asymmetry of the reaction. This 
happens because of a decrease in the participant zone. 
Although the role of the Coulomb interaction decreases 
with energy, its effect, however, remains constant ( 20%) 
with the asymmetry of the reaction. Because of the 
presence of the Coulomb forces, the nuclear matter 
breaks into smaller pieces/free nucleons. 



FIG. 3: (Color online) The mass asymmetry r\ variation of 
the production of A max , free nucleons, LMF's and IMF's with 
and without the Coulomb effect at two different energies at 
Ec. m . = 50 MeV/nucleon and E c . m . = 250 MeV/nucleon and 
impact parameter at b = 0.3. 

To understand this aspect further, we display in Fig. 4, 
the saturation density of the reaction obtained at 
200 fm/c. This density remains quite high for larger 
asymmetries. This is caused by the least amount of 
destruction of nuclear matter at larger asymmetries, 
which leads to higher nucleonic densities. As a result, 
one sees a heavier A max compared to the smaller 77 
values. Because of the heavier A max , the emission of the 
nucleons is also smaller. 

In Fig. 5, the charge distribution is displayed as a func- 
tion of the fragment charge using S (upper panel) and H 
EOS (lower panel). The symmetric reactions (?7=0) lead 
to enhanced emission of nucleons and LMFs compared 
to asymmetric reactions, where incomplete fusion events 
or deep inelastic events are dominant. We can also say 
that a large asymmetric reaction leads to few nucleonic- 
transfer processes. One also notices that the slope of the 
distribution becomes steeper with H EOS compared to S 
EOS. Because of the enhanced binary collisions between 
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FIG. 4: (Color online) The variation of density with mass 
asymmetry at an impact parameter b — 0.3. The different 
symbols are at E c . m . ~ 50 MeV/nucleon and E c . m . — 250 
MeV/nucleon with S and H EOS. 



the nucleons for nearly symmetric nuclei, the emission of 
fragments is suppressed. Clear systematics can be seen 
in the production of fragments with the asymmetry of 
the reaction. 

In Fig. 6, the variation of the multiplicity of LMFs is dis- 
played as a function of the center-of-mass energy for vari- 
ous asymmetric reactions using H and S EOS. Because of 
more compression, the nearly symmetric reaction drives 
matter into the participant zone and, as a result, more 
lighter fragments are emitted. 

In Fig. 7, the variation of the multiplicity of IMFs 
is displayed as a function of the center-of-mass energy 
E c . m .. This happens because of the fact that the system 
suffers less compression; and, hence, less numbers of 
IMFs are produced. One notices several interesting 
points: The nearly symmetric collision leads to a well- 
known trend (i.e., the maximum emission occurs around 
100 MeV/nucleon). This trend, however is not shown 
by the large asymmetric reactions where we donot see 
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ent asymmetries between r\ — 0.2 to 0.7 at E c . m . = 150 
MeV/nucleon and impact parameter at b = 0.3. The upper 
and lower panels are shown with S and H EOS, respectively. 



any sharp rise or fall; and, furthermore, a flat plateau 
is obtained at much higher incident energies compared 
to nearly symmetric nuclei. Therefore, experiments are 
needed to verify this prediction. 



IV. CONCLUSION 

A systematically theoretical study is presented for 
the asymmetric colliding nuclei, which use a variety of 
reactions that employ different EOS as well as incident 
energies. We envision an interesting outcome for large 
asymmetric colliding nuclei, although nearly symmetric 
nuclei depict a well-known trend of rising and falling 
with peak around E= 100 MeV/nucleon. This trend, 
however, is completely missing for large asymmetric 
nuclei. In conclusion, experiments are needed to verify 
this prediction. 
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